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Abstract 

We prove that solutions of stochastic differential equations driven 
by fractional Brownian motion for H > 1/2 define flows of homeomor- 
phisms on R''. 

AMS Subject Classification: 60H05, 60H07 

1 Introduction 

Suppose that = {B^ ,t > 0} is an m-dimensional fractional Brownian 
motion with Hurst parameter H G (0, 1), defined in a complete probability 
space {Q,J-,P). That is, the components B^'^ , j = 1, . . . ,m, are indepen- 
dent zero mean Gaussian processes with the covariance function 

1 ''.2H , 2s u „|2H 



RHit,s) = ^[t'^+s'^-\t-sr). (1.1) 

For if = i, the process B is an m-dimensional ordinary Brownian motion. 
On the other hand, from ()1.1() it follows that 



EiBf'^-B^^^ ) = \t-s 

As a consequence, the processes B^'^ have stationary increments, and for 
any a < H we can select versions with Holder continuous trajectories of 
order a on a compact interval [0, T]. 

This process was first studied by Kolmogorov in [7j and later by Man- 
delbrot and Van Ness in ^T], where a stochastic integral representation in 
terms of an ordinary Brownian motion was established. 

The fractional Brownian motion has the following self-similar property: 
For any constant a > 0, the processes \^a~^ B^^,t > 0} and {^B^ ,t > 0} 
have the same distribution. For H = \ the process B^ has independent 
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increments, but for H ^ this property is no longer true. In particular, 
\i H > ^, the fractional Brownian motion has the long range dependence 
property, which means that for each j = 1, . . . , m 

oo 
n=l 

The self-similar and long range dependence properties make the frac- 
tional Brownian motion a convenient model for some input noises in a vari- 
ety of topics from finance to telecommunication networks, where the Markov 
property is not required. This fact has motivated the recent development of 
the stochastic calculus with respect to the fractional Brownian motion. We 
refer to for a survey of this topic. 

In this paper we are interested in stochastic differential equations on 
driven by a multi-dimensional fractional Brownian motion with Hurst 
parameter > ^, that is, equations of the form 

^ j-t pt 
Xl = Xi + y^ a^'^{s,X,)dB^'^ + / b\s,Xs)ds, (1.2) 
Jo Jo 

i = l,...,d. The stochastic integral appearing in H1.2() is a path- wise 
Riemann-Stieltjes integral. In fact, under suitable conditions on a, the 
processes (t{s,Xs) and Bf have trajectories which are Holder continuous 
of order strictly larger than ^ and we can use the approach introduced by 
Young in 17 . A first result on the existence and uniqueness of a solution 
for this kind of equations was obtained by Lyons in 8 , using the notion 
of p-variation. On the other hand, the theory of rough path analysis intro- 
duced by Lyons in |Sj (see also the monograph by Lyons and Qian ^jOj ) , has 
allowed Coutin and Qian |3j to establish the existence of strong solutions 
and a Wong-Zakai type approximation limit for the stochastic differential 
equations of the form (|1.2|) driven by a fractional Brownian motion with 
parameter H > ^. In i3| sufficient conditions on the vector fields b and a 
are given to ensure existence and uniqueness of the solution of (|1.2|1 even 
when vector fields do not commute. 

In jl8j Zahle has introduced a generalized Stieltjes integral using the 
techniques of fractional calculus. This integral is expressed in terms of frac- 
tional derivative operators and it coincides with the Riemann-Stieltjes in- 
tegral fdg, when the functions / and g are Holder continuous of orders 
A and fi, respectively and A + > 1. Using this formula for the Riemann- 
Stieltjes integral, Nualart and Ra§canu have obtained in jTl] the existence of 
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a unique solution for the stochastic differential equations ()1.2() under general 
conditions on the coefficients. 

Later on, Nualart and Saussereau ^] have studied the regularity in the 
sense of Malliavin Calculus of the solution of Equation (|1.2p , and they have 
established the absolute continuity of the law of the random variable Xt 
under some non-degeneracy conditions on the coefficient a. 

The main result of this paper is the flow and homeomorphic properties of 
X as a function of the initial condition x. Since the solution of ()1.2() is defined 
path- wise, ordinary (i.e., deterministic) methods are in use here. Namely, we 
use the estimates found in 14 and approximate fractional Brownian motion 
by a sequence of regular processes to prove that the solution {Xj.t{x), < 
r <t <T, X eW^} oi 

Xrt{x)=x+ I a{s,Xrs{x))dB"{s)+ I b{s,Xrsix))ds, 

defines a flow of R'^-homeomorphisms. 

The paper is organized as follows. Section 2 contains some preliminaries 
on fractional calculus. In Section 3 we review some results on the properties 
of the solution of Equation H1.2|) and we establish some continuity estimates 
as a function of the initial condition and the driven input, which are needed 
later, finally in Section 4 we show that Equation 1)1.2(1 defines a flow of 
homeomorphisms. 



2 Preliminaries 

Let a,b € , a < b. Let f £ L} (a, b) and q > 0. The left-sided and right-sided 
fractional Riemann-Liouville integrals of / of order a are defined for almost 
all X G (a, b) by 

(x) = -1- r (x - yf-' f (y) dy 
r (a) J a 

and 

r(a) Jx 

respectively, where (—1)"" = e~*'^" and T (a) = r°'~^e~^dr is the Euler 
function. Let I"_^_{LP) (resp. I^_{LP)) the image of U'{a,b) by the operator 
/"+ (resp. I^J). U f £ 12+ (LP) (resp. / € (LP)) and < a < 1 then 
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the Weyl derivative 

^.\/ (-) - + ^^^dy) (2.3) 



(2.4) 

is defined for almost all x G (a, 6) . 

For any < A < 1, denote by C^(0,T) the space of A-H61der continuous 
functions / : [0,T] — > M, equipped with the norm \\f\\^ + ||/||;^, where 

sup ^m^. 

0<s<t<T [t — S) 

Recall from jl6j that we have: 

• If a < - and g = ^ then (LP) = 7° (LP) C L« (a, b) . 

p 1 — ap 

• If a > ^ then/^+ {LP) U (L^) C C^'p (a, 6) . 

The linear spaces (L^) are Banach spaces with respect to the norms 

ll/ll/^^(LP) = II/IIlp + ll^a+/|Lp ~ ||^a+/|Lp > 

and the same is true for (LP). 

Suppose that / G C^{a,b) and g £ C'^{a,b) with A + /i > 1. Then, 
from the classical paper by Young 17 , the Riemann-Stieltjes integral fdg 
exists. The following proposition can be regarded as a fractional integration 
by parts formula, and provides an explicit expression for the integral fdg 
in terms of fractional derivatives (see '18^^ ) . 

Proposition 2.1 Suppose that f G C^{a,b) and g G C'^{a,b) withX+jj > 1. 
Let X > a and /x > 1 — a. Then the Riemann Stieltjes integral jj* fdg exists 
and it can be expressed as 

f fdg = (-1)" /' D^^f (t) Dlz'^gb- (t) dt, (2.5) 

J a J a 

where gi,^ (t) = g{t) -g{b). 
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In ^Sj Zahle introduced a generalized Stieltjes integral of / with respect 
to g defined by the right-hand side of (|2.5() . assuming that / and g are 
functions such that gib—) exists, / G I^^ (L^) and gi,^ E -^6-" i^*^) 
some p,q > 1, 1/p + 1/q < 1, < a < 1. 

Let a < i and d G N*. Denote by Wo"'°°(0, T; M'^) the space of measur- 
able functions f : [0,T] —>■ such that 



l/ILoo- -p (I/WI+/ .I7' ds]<^. 



We have, for all < e < a 

C7"+"(0,T;M'^) C Wo"'°°(0,r;M'^) C C7"-^(0, T; R'^). 
Denote by W^""'°°(0, T; R™) the space of measurable functions g : [0, T] 



r>m 



such that 



\git)-g{s)\ , 



Clearly, for all e > we have 

Moreover, if belongs to W^~°'°°(0, T; M™), its restriction to (0,t) belongs 
to //r"(L°°(0,t;M")) for ah t and 

^"(^)-=wT^ ^ \i^l-"9t-){s)\ 

r(l - a) o<s<t<T 

^r(i-l)r(a) Nli-.oo,T < 00. (2.7) 

The integral Jq /d(7 can be defined for all t G [0, T] iig belongs to W^~"'°°(0, T) 
and / satisfies 

Furthermore the following estimate holds 



fdg 







<A.(5)ll/IL,i 
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3 Stochastic differential equations driven by a fBm 



We are going to consider first the case of a deterministic equation. Let 
< a < ^ be fixed. Let g G l^y~"'°°(0, T; M™'). Consider the deterministic 
differential equation on 



^0 Jo 



l,...,d, where xq G M*^, and the coefficients a^'^,¥ : [0,T] x 



measurable functions. Set a 



' dxm 



are 



, , and for a matrix 

2 



dxm 
2 



and a vector y = {u j^y^i denote |^| = X^ijn*'"'! 

\y\^ = J2i\y , 

Let us consider the following assumptions on the coefficients. 

(HI) a{t,x) is differentiable in x, and there exist some constants < 
P,6 < 1, Ml, M2, M3 > such that the following properties hold: 



i) Lipschitz continuity 

\a{t,x) - a{t,y)\ < Mi\x - y\, VxGM'^, VtG [0, T] 

ii) Holder continuity 

\d^^ait,x) -d^^a{t,y)\ < M2\x - y\^ , 

V|x| , \y\ £ M"^, Vt G [0,r] ,i = 1,... ,d, 
in) Holder continuity in time 

\a{t,x) - a{s,x)\ + \dx^(7{t,x) - dx^(T{s,x)\ < Ms]* - sf 

Vx G M'^, Vt,s G [0,T] . 



(H2) There exists constants Li,L2 > such that the following properties 
hold: 

i) Local Lipschitz continuity 

\b{t,x) -blt,y)\ <Li\x-y\, V|x|,|y| G M"^, Vi G [0,r], 

Linear growth 



Set 



\b{t,x)\ < L2(l + |a;I), Vx G R'^, Vt G [0,r] 



ao = mm < -, p, ■ 



2 ' 1 + 5 , 

The following existence and uniqueness result has been proved in |14j . 
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Theorem 3.1 Suppose that the coefficients a{t,x) and b{t,x) satisfy as- 
sumptions (HI) and (H2). Then, if a < uq there exists a unique solution of 
Equation in the space C^~°' (O, T; R'^) . 

Actually, these conditions can be slightly relaxed. For instance, the 
Holder continuity of the partial derivatives of a and the Holder continuity 
of the coefficient b may hold only locally (see |2] for the details). 

We now state two theorems which are consequences of the estimates 
found in ^1]. 

For any A > we introduce the equivalent norm in the space Wq'°^ (O, T; M.'^^ 
defined by 

ll/lkA= -P (m\+f '^f "{ff' ds 

tG[0,T] V Jo [i- S> 

Theorem 3.2 Let us denote £,t{xo) the solution of at time t with initial 
condition xq. Fix R > 1. Then there exists a constant C such that for any 
xq and xi in the hall B{0, R) = {x £ W^, \x\ < R} and for any A > 



Rexp(C{l + Aaig))^ 



we have 



Uixo) - C{xi)\\a,x< (1 - ^exp {^C (1 + A„((7))^ j A^°-^j [xq - xi\. 

Proof. It is proved in that there exists a constant Ci such that 

1 

if Aq = Ci (1 + Aq,(5()) i-2c« then for any initial condition xq in the ball of 
radius R we have 

||^(xo)IU,Ao<2(l + l^o|)<4i?. (3.9) 
Given a function / G Wq''^{0, T; M'^) we define as in jl4j 

(/)= [\{s, f{s))ds, 
Jo 

C''Hg.f)=j'<T(s. HB))dg(s), 



rG[0,T] Jo 
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If /, /i e Wq'°°{0, T; R'^) (see [H]) there exist constants C2 and C3 such that 



F^^\f) - F^^\h) < CaA"'^ 11/ - /i 

a, A 

G(-) (5, /) - G(-) (5, /i) ^ < C3A.(<7)A2"-i 

a, A 



la,A ' 



(3.10) 



x{l + A{f) + A{h))\\f-h\U, (3.11) 



^('^^5,/) , <C4A„(5)A2--i(i + ||/| 



a, A 



cj,A 



(3.12) 



Also, if / € Wq ' (Oi ^! ^ ) /i is a bounded measurable function, then 



00/ ' 



F^'Hh) <C5(1 + 

1 — Q 

G^'^Uff,/) <C6A,(5)(l + 
We have the following estimate 



(3.13) 
(3.14) 



^{C{xo)) = sup 

re[0,T] ./O 



(r — s) 



a+l 



ds 



rp5-a(l+S) 

— -^iie(xo)iii_„, 



5-a(l + (5) 
and using ()3.13|) . ()3.14|) . and ()3.9() we obtain 



(3.15) 



+ 



1-0 



G('^) (e(xo)) 



< |xo| + C5 (1 + llC(xo)ll^) + A,(<7)C6 ( 1 + |!C(xo)||, 
<2C7e^«^i?(l + A„((7)) 

= 2C7 exp (tCi (1 + A,((7))i^ ) R{1 + A„(5)). 



Hence, from (|3.15|) we get 

A(e(xo)) < Cs exp frCi (1 + Aa{g))T^ ) R{1 + A,(5)). 



(3.16) 



Using ()3.1U() . (|3.11() . and 1)3. 16() we obtain for xq and xi in the ball of radius 
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R 



+ 



G('^^(5,e(xo))-GM(5,e(xi)) 



< \xo-xi\ + CiX''-^U{xo)-axi)L,x 

+ C72A„(5)A2"-i (1 + a (C(xo)) + A (e(xi))) ll^(xo) - 

< |xo-xi[ + CiA°-i||C(xo)-e(xi)L,A 

+ C72A,(5)a2"-i (l + 2Cs exp (tCi (1 + Aa{g))^) R{1 + A^ig)) 
X ||^(xo)-e(2;i)L,A- 

As a consequence, 

ll?(xo)-e(xi)L,A < |xo-xi|+A2-i (exp (C7(l + A„(5))^) i?) Uixo) - axi)L,x 
for some constant C, which imphes the result. ■ 
Theorem 3.3 The map 

where ^ is i/ie solution of it A'. uiii/i xq € ^(0, i?), is continuous. Namely, 
for A > large enough we have 

Proof. With the previous notations, we have 

and 

^(/i) = xo+FW(?(/i))+GM(/i, m)- 

Thus, 

e(<7)-ew = F('')(e(<7))-F(')(^W) 

+ {g -h,a9)) + (/^, - (/i, . 
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According to H3.12() . for A sufficiently large, 

+ C4A2"-l(l+A,(5-/l))||e(5)|la,A 

+ CsX^''~'A^{hm9)-ah)\\a,x- 

Hence the result. ■ 

Consider Equation 1)1.2^ on M*^, for t € [0, T], where Xq is a d-dimensional 
random variable, and the coefficients a^'^,b^ : Q x [0,T] x M'^^M are mea- 
surable functions. 

Theorem 3.4 Suppose that Xq is an W^-valued random variable, the co- 
efficients a{t,x) and b{t,x) satisfy assumptions (HI) and (H2), where the 
constants might depend on to, with (3 > 1 — H, 6 > 1/H — 1. Then if 
a G (1 — H,ao), then there exists a unique stochastic process 

X G L° (n, P ; VFo"'°°(0, T; M^)) 

solution of the stochastic equation hl.l3\) and, moreover, for ¥-almost all 
w G J7 

X {lo, •) = (X* (a;, G C''^ (o, T; M'^) . 

Consider the particular case where 6 = and a is time independent, that 
is, ^ 

Xt = Xo+ [ a{Xs)dBf. (3.17) 
Jo 

By the above theorem this equation has a unique solution provided a is 
continuously differentiable, and a' is bounded and Holder continuous of order 
6>jj-l. 

In Hu and Nualart have established the following estimate. Choose 
9 G Then, the solution to Equation ()3.17|) satisfies 

sup \Xt\ < 21+'=^(II'^'II-^I-(0)I)I|b«I|J'"'(|j^q| + i) . (g.ig) 

0<t<T 

If a is bounded and ||cj'|| ^ this estimate can be improved in the following 

way 

sup \Xt\<\Xo\ + k\\a\\oo(T^\\B"\\l "^na'lffllB^lllY (3.19) 

0<t<T \ J 
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These estimates improve those obtained by Nualart and Ra§canu in 
based on a suitable version of Gronwall's lemma. The estimates ()3.18|) and 
H3.19() lead to the following integrability properties for the solution of Equa- 
tion (1X171) . 

Theorem 3.5 Consider the stochastic differential equation \'J.n\j , and as- 
sume that E{\Xq\p) < oo for all p > 2. If a' is bounded and Holder contin- 
uous of order 6 > — 1, then 

E \ sup \Xt\n < oo (3.20) 
yo<t<T J 

for all p > 2. If furthermore a is bounded and E (exp(A|Xo|''')) < oo for any 
A > and 7 < 2H , then 

E ^expA ^ sup \Xt[^ < 00 (3.21) 

for any A > and 7 < 2H . 

In l15_ Nualart and Saussereau have proved that the random variable Xt 
belongs locally to the space 15°° if the function a is infinitely differentiable 
and bounded together with all its partial derivatives. As a consequence, 
they have derived the absolute continuity of the law of Xt for any t > 
assuming that the initial condition is constant and the vector space spanned 
by {(a*'(xo))i<i<d, 1 < J < m} is M'^. 

Applying Theorem 13.51 Hu and Nualart have proved in [B| that if the 
function a is infinitely differentiable and bounded together with all its partial 
derivatives, then for any t G [0, T] the random variable Xt belongs to the 
space B°°. As a consequence, if the matrix a(x) = aa'^{x) uniformly elliptic, 
then, for any t > the probability law of Xt has an C°° density. In a recent 
paper, Baudoin and Coutin P have extended this result and derived the 
regularity of the density under Hormander hypoellipticity conditions. 

4 Flow of homeomorphisms 

Let vr = {0 = to < ^1 < ■ ■ • < = T} be the uniform partition of the interval 
[0,r]. That is tfc = ^, /c = 0, ...,n. We denote by S"-^ the polygonal 
approximation of the fractional Brownian motion defined by 

n-l 
k=0 
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In order to get a precise rate for these approximations we will make 
use of the following exact modulus of continuity of the fractional Brownian 
motion. There exists a random variable G such that almost surely for any 
s, t G [0, T] we have 



-B^ \ < G\t 



HI 



\/log(|t 



(4.22) 



Fix 9 < H. We have the following result, which provides the rate of conver- 
gence of these approximations in Holder norm. 



Lemma 4.1 There exist a random variable Ct,i3 such that 



C»(0,T;R'") 



(4.23) 



Proof. To simplify the notation we will assume that m = 1. Fix < 
s < t < T and assume that s £ [ti,ti^i] and t E Let us first 

estimate 

hi{s,t) 



1 



it - sY 



IB-''' - B^ 



[Bf^'-B^) 



If t — s > ^ , then using (|4.22j) we obtain 

\hi{s, t)\ < T-^n^ \ Bll -B? + ^{t- t,) - B^ 



+ 



< 4Gr-^+^n-^+Vlog {n/T). 

If t — s < ^, then there are two cases. Suppose first that s, t G [ffc, ifc+i]- In 
this case, if n is large enough we obtain using ()4.22|) 

\Bi^-B^\ n\Bi' -Bf^l 

II,/ J.\\ ^ ' t , " I *fc ' /j. \ 

\hi{s,t)\ < — + - — — {t - s) 



{t - sf 



< G\t - s|^-Vlog|i-sh^ + Gr-i+^0og(n/r) n^-'^it - sf-*^ 

< 2Gr-^+^n-^+Vlog (n/T). 



On the other hand, if s G [tk-i,tk] and t G [tk,tk+i] we have, again if n is 
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large enough 



\hi{s,t)\ < 



it - s)' 

{<-i?f-^(..-.)(<-<_^)} 



< 



1 



< 



it - sY 

G 



\Bf'-B^\ + ^it-s){\Bll-Bflj + \B!^ 



k + l 



B 



\t-sf^\og\t-s\-^ + 2{t-s) (5)" ^/^og{n/T) 



H 



< 3Gr-^+^n-^+Vlog (n/T). 
This proves ■ 

Corollary 4.1 For any a £ {1 — H, 1/2), we /ia?;e; 



supAa(S"'^) < +00 and hm A„(B"'^ _ S«) = 0. 

n— >+oo 

Proof. Choose r] > in such a way that 1 — a + t] < H . According to 
(fTK|l and (fT7|l . we have 

Aa(-B"''"^) < CfJi?"''^||(7i-a+r,(o,T;R'") 

and 

AaiB"-'^ - B^) < CriWB"-'^ - B^\\cl-a+v(o,T;K'^)- 

Then, Lemma 14.11 imphes that the sequence B"'' converges to B^ in the 
norm of C^~"~^^{0,T;MJ^) which yields the results. ■ 

Consider for any < r < t < T and any natural number n > 1 the 
following equations 



X^t{x) = x+ I a{s,Xj;,{x))dB'''^{s)+ I b{s, Xl^,{x))ds, (4.24) 



and 



Y;i{x)=x+ I a{s,Y:i{x))dB^'''{s)+ I h{s X^t{^))ds . (4.25) 



We know from standard results on ordinary differential equations that for 
any n > 1, 



1. Equations (|4.24|) and (|4.25|) have a unique solution. 
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2. For any x G M'', for any < r < r < t < T, X'!^t{X^^{x)) = X^^t^x). 

3. For any x G M'^, for any < r < r < t < T, YJ'^iY^Kx)) = Y;1{x). 

4. The maps (x i— > X^^{x)) and {x i-^ Y^{x)) are M'^-homeomorphisms 
inverse of each other: 

= X and (x)) = x. 

We are then in position to prove our main theorem: 

Theorem 4.1 Assume that Hypothesis (HI) and (H2) hold. Then, claims 
^\^li^and^ also hold for the equations 

Xrt{x)=x+ j a{s,Xrs{x))dB^{s)+ j b{s,Xrs{x))ds, (4.26) 

J r J r 

and 

Yrt{x)=x+ [ a{s,Y,t{x))dB^{s)+ [ b{s,Yst{x))ds. (4.27) 

Jr J r 

Proof. Point ^ is proved in ^1] . As to the second claim, proceed as 
follow: 

X,",(X,",(X)) - XrtiXrrix)) = X!^,{Xl\{x)) - X!^,{Xrr{x)) 

+ iX!^t-Xrt)iXrrix)). 

Fix e > and a such that 1 — H < a < ^. Fix a trajectory a; G il. 
Choose no so that Ka{B^'^ — B^) < e for all n > uq and choose A such 
that A^°~^C2 sup„ Aa(-B"'''^) < |. Then, according to Theorem 13.31 for any 

n > no, 

\ 2Q! — 1 II II 

\\X,. - Xr.\\a,X < ^_c^^2a-l (1+A„(ijn,/^)) ^ ^^1^^ l|A..|U,Ae. 

Hence, for n > uq, 

\{X!^t - Xrt){XrAx))\ < ce. 

The convergence of X" implies that there exists R such that for any r G [r, t] 
and for any n > no, X^^{x) G B{0,R). Then, Theorem 13.21 implies that for 
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A large enough 

\X^,{X:^,{X)) - X^,{Xrr{x))\ 

< f 1 - iiexp ( C ( 1 + supA„(5 



1 \ X -1 

n,H\ \ 1 \2q-1 I 



X |X;,(x)-X,,(x)| 

< C\XI^^{X) - Xrr(x)\. 

We have thus proved that 

0= hm X:^t{Xl^,{x)) - XrtiXrr{x)) 

= hm X^,{X) - Xrt{Xrr{x)) 

n— >+oo 

= Xrtix) - Xrt{Xrr{x)). 

Other points are handled similarly. ■ 

Acknowledgement 1 This work was carried out during a stay of Laurent 
Decreusefond at Kansas University, Lawrence KS. He would like to thank 
KU for warm hospitality and generous support. 



References 

[1] F. Baudoin and M. Hairer: A version of Hormander's theorem for the 
fractional Brownian motion. Preprint, 2006. 

[2] L. Coutin, P. Friz, and N. Victoir: Good rough path sequences 
and applications to anticipating and fractional stochastic calculus. 
|arXiv:m ath.PR/0 501197| 

[3] L. Coutin and Z. Qian: Stochastic analysis, rough path analysis and 
fractional Brownian motions. Probab. Theory Related Fields 122 (2002) 
108-140. 

[4] L. Decreusefond and A.S. Ustiinel: Stochastic analysis of the fractional 
Brownian motion. Potential Anal. 10 (1999) 177-214. 

[5] D. Feyel and A. de La Pradelle: On fractional Brownian processes. Po- 
tential Anal. 10 (1999) 273-288. 



15 



[6] Y. Hu and D. Nualaxt: Differential equations driven by Holder continu- 
ous functions of order greater than 1/2. Preprint, 2006. 

[7] A. N. Kolmogorov: Wienershe Spiralen und einige andere interessante 
Kurven im Hilbertschen Raum. C. R. (Doklady) Acd. Sci. URSS (N. S.) 
26 (1940) 115-118. 

[8] T. J. Lyons: Differential equations driven by rough signals. I. An exten- 
sion of an inequality of L. C. Young. Math. Res. Lett. 1 (1994) 451-464. 

[9] T. J. Lyons: Differential equations driven by rough signals. Rev. Mat. 
Iberoamericana 14 (1998) 215-310. 

[10] T. J. Lyons and Z. Qian: System control and rough paths. Oxford 
Mathematical Monography, (2002). 

[11] B. B. Mandelbrot and J. W. Van Ness: Fractional Brownian motions, 
fractional noises and applications. SIAM Review 10 (1968) 422-437. 

[12] D. Nualart: Stochastic integration with respect to fractional Brownian 
motion and applications. Contemp. Math. 336 (2003) 3-39.. 

[13] D. Nualart and E. Pardoux: Stochastic calculus with anticipative inte- 
grals. Probability Theory and Related Fields 78 (1988) 535-582. 

[14] D. Nualart and A. Rascanu: Differential equations driven by fractional 
Brownian motion. Collect. Math. 53 (2002) 55-81. 

[15] D. Nualart and B. Saussereau. Malliavin calculus for stochastic differ- 
ential equations driven by a fractional Brownian motion. Preprint, 2005. 

[16] S.G. Samko, A. A. Kilbas, and O.I. Marichev: Fractional integrals and 
derivatives, Gordon and Breach Science, 1993. 

[17] L. C. Young: An Inequality of the Holder Type Connected with Stieltjes 
Integration, Acta Math. 67 (1936) 251-282. 

[18] M. Zahlc: Integration with respect to fractal functions and stochastic 
calculus, I. Probability Theory and Related Fields 111 (1998) 333-374. 



16 



